Magnetic resonance spectroscopy and characterization of magnetic phases for spinor 

Bose-Einstein condensates 
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The response of spinor Bose-Einstein condensates to dynamical modulation of magnetic fields is 
discussed with the linear response theory. As an experimentally measurable quantity, the energy 
absorption rate (EAR) is considered, and the response function is found to access quadratic spin 
correlations which come from the perturbation of the quadratic Zeeman term. By applying our 
formalism to spin-1 condensates, we demonstrate that the EAR spectrum as a function of the 
modulation frequency is able to characterize the different magnetically ordered phases. 
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Introduction. — Ultracold Bose atoms with spin de- 
grees of freedom have been attracting interests as a 
class of quantum fluids accompanying non-trivial spin or- 
ders and topological spin textures, in contrast with spin- 
less Bose-Einstein condensates (BECs). In BECs with 
spins, so-called spinor BECs, spin rotational symmetry 
allows spin dependent interactions, and the number of 
the independent interactions increases with spin degrees 
of freedom of atoms. Furthermore, such spin dependent 
interactions cause various ground states. However, in 
addition to exploring the properties of those nontrivial 
spin orders, it is also important to specify them exper- 
imentally. Thus the development of the measurement 
techniques to capture complicated orders is a challenge 
for the study on spinor BECs. 

There exist the several experimental techniques for 
spinor BECs: the spin texture imaging by the sep- 
aration of the different spin components by magnetic 
field gradient, the spatially- resolved magnetometer by 
observing Larmor precession 04^, and measurement of 
phase transitions and spin dynamics by dynamical con- 
trol of quadratic Zeeman (QZ) shift [7hTo|. In addition, 
measurement of birefringence by the off-resonant light 
scattering technique has been analysed, which can ac- 
cess structure of spin textures such as nematicity and 
Skyrmions. [1, [ll| These probes are mainly used to cap- 
ture static properties of systems, and it is not easy to 
capture excitation spectra involving spin degrees of free- 
dom. 

For spinor BECs in the presence of uniform magnetic 
fields, the QZ shift in addition to the linear Zeeman (LZ) 
shift emerges due to hyperfine couplings between a nu- 
clear and an electron spin. The LZ shift in experi- 
ments is considered to be negligible because of spin con- 
servation, and indeed the magnetization of spinor BECs 
is preserved at least within the limit of accuracy of experi- 
mental errors, [l^l Thus the most relevant effect induced 
by the magnetic field is the QZ shift, and furthermore 
what is important is that this QZ coupling is an experi- 
mentally controllable parameter. [gI. [sIIqI [l^ . 

In this paper, motivated by such possible control of 
the QZ coupling, we consider magnetic resonant spectra 



as a response to dynamically modulated magnetic fields. 
As an experimentally measurable quantity, we focus on 
energy absorption rate (EAR), and formulate this EAR 
spectrum as a function of the modulation frequency of 
magnetic fields by using linear response theory. In ad- 
dition, to demonstrate the potential application of this 
spectroscopy, we analyse the spectra for uniform spin- 
1 BECs with Bogoliubov theory [l^ [l^, and the con- 
sequent spectra are found to exhibit different behaviors 
in each phase. Finally we also consider the cases in the 
presence of trap potential, and conclude that the ordered 
phases of the trapped systems remains distinguishable 
from the low-frequency behavior of the EAR spectrum. 

Formalism. — We start with general spin-i^ Bose atom 
systems under an uniform magnetic field. Let us sup- 
pose the many-body static Hamiltonian including Zee- 
man couplings to be Hq. In this paper, we restrict our- 
selves that the magnetic field is applied along z axis, and 
the Hamiltonian is invariant under spin rotation around 
z axis, which is a general setup of the experiments. 

In the presence of the dynamically modulated mag- 
netic filed such as h -\- 6hcos{ujt/2), the system should 
be described by the time-dependent Hamiltonian H{t) = 
Hq + V{t), and the perturbation is represented as 

V{t) = {Sp) cos(a;i/2) -Fl + (Sq) cos^ {ujt/ 2) Tq, (1) 

where the first and second terms mean modulation of 
the LZ and QZ couplings, respectively. The coupling 
constants, Sp and Sq, are proportional to Sh and {5hY , 
respectively. The LZ and QZ operators are represented 
as 



(2) 
(3) 



where ^I" — [ippjipp-i, ■ ■ ■ ^ip-p)'^ denotes a spinor bo- 
son field, and F = (F^, F^ , F^) is a spin-F matrix. 

In experiments, the EAR can be measured through 
the time of flight image. Assuming the energy scale 
of the periodically modulated perturbation to be small 
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enough [17|, the dynamics is well-described with the 
linear response theory. Then, the EAR is defined as 



^{'^) = 2i7iJ It rft ' where (• • •) denotes the 
statistical average by H{t). Thus the EAR is derived as 



Ricu) 



2h 



ujlm[x (uj) 



(4) 



where x^{uj) is the Fourier transformed retarded cor- 
relation function of the perturbation ([T]): x^(i) — 
—i {\y [t) ,V {Q)]) Q averaged by H^. Since is assumed 
to possess the spin rotational symmetry around z axis, 
[JL,-ffo] = [-^Lj-^q] = 0, and thus the retarded correla- 
tion function is reduced to 
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FIG. 1. (Color online) The MF phase diagram of spin-1 BECs. 
The phase boundary between the BA and LP phase is given 
by q = 2n|ci|. The thick-solid and dashed boundaries denote 
a first and second order phase transition, respectively. 



:[-^qW,-^q(o)] 



(5) 



Namely, the system is insensitive to the dynamic modu- 
lation of the LZ coupling. 

EAR for spin-1 BECs. — Let us demonstrate the EAR 
spectrum Q to allow for characterizing spin-ordered 
phases. We consider spin-1 interacting bosons without a 
trap, which undergo BEC in the low-temperature regime. 
The Hamiltonian to be considered is given by 



-ffn- dr 



, (6) 



where 4' = {ipi,ipo,ip-i)'^ , M denotes a mass of atoms, 
and Co and ci mean the density and the spin exchange 
interaction, respectively. For ^^Na and *^Rb atoms, the 
coupling ci is taken to be positive and negative value, 
respectively. Let us impose ncg ^ ''T-|ci|, \q\, where n is a 
atom density to meet the experimental conditions. Since 
the EAR is independent of the LZ coupling modulation 
as discussed above, instead of (HI we here suppose the 
modulation perturbation of the magnetic field to be 



V{t) = [Sq) cos(L^t)J-Q. 



(7) 



The mean-field (MF) analysis, in which the field ^ 
is replaced by a MF spinor order parameter ^ optimiz- 
ing the Hamiltonian (jH]), leads to the following ground 
states [13 as shown in Fig.[T] (i) the ferromagnetic (FM) 
phase ^FM = (1,0,0)-^ for ci,q < 0, (ii) the longitudinal 
polar (LP) phase ^lp = (0,1,0)"'" for ci,q > or for 
ci < and q > 2n\ci\, (iii) the transverse polar (TP) 
phase Ctp (1/^2,0, l/V2)'^ for ci > and g < 0, and 
(iv) the BA phase |ba = (sin 6'/%/2, cos 6*, sin6'/%/2)^, 
where sin6' = ^(1 — q)/2 with q = q/{2n\ci\), for ci < 
and < g < 2n|ci|. 

The correlation function x^{uj) is calculated with the 
Bogoliubov theory, by the replacement ^ = V^^a + 
(a = FM, LP, TP, BA) where Nq is the condensate atom 



number. Then the Bogoliubov Hamiltonian is repre- 
sented as Hq « Emf + Hcs with the MF energy Emf 
and 

HcS = ^ 4 (efc + nciF ■ {F -F) + q{{F' f - W?)) 
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(^Sl -Sk + Sk- S-k + H.c 



(8) 



where dk = Jdre^''"^S^ = (aifc, aofc, a_ifc)-^ is a Fourier 
transform of the fluctuation of the spinor. Dk — £,^dk 
and Sk — Fdk mean a density and spin fluctuations, 
respectively, and X — ^^X£^ for a spin operator X de- 
notes the MF average of spin matrices. In this represen- 
tation, the QZ operator is expressed as 



(9) 



FM phase. — The MF spinor leads to the following 
density and spin fluctuation as Dk = — aik, — 
o.ok/V2, and S^. — iaok/V2. The eff'ective Hamilto- 
nian diagonalized by the Bogoliubov transformation is 
given [16] as 



i?™(fc)4(fc)dF(fc) + Ei^{k)fi{k)f,xk) 



Elf{k)fl^{k)f,,y{k) 



(10) 



where E™{k) = J ek[ek + 2n{c^ - \c,\)], E^^{k) = 
Ck + 2|ci|, and E™{k) ^ Ck + \q\. According to the 
Bogoliubov transform, the QZ operator in the FM phase 
is written as 



J-^^ = Const. 



k^O 



fljk)Uy{k). 



(11) 



The perturbation ^™ commutes with the Hamilto- 
nian (1101) . and immediately the response function x^(cli) 



3 




FIG. 2. (Color online) The EAR spectra as a function of modulation frequency for the different QZ couplings q = q/2n\ci\: (a) 
the LP phase for ci, g > or for ci < and q > 2n|ci|, (b) the TP phase for ci > and q < 0, and (c) the BA phase for ci < 
and < q < 2n\ci\. The spectra (a) and (b) show qualitatively the same behavior, but the intensity differs quantitatively due 
to the different number of the accessible gapful-spin modes. The spectra (c) apart from q — 1 show the abrupt enhancement 
of the EAR around hcj ~ 2n|ci|, which indicates the appearance of the contribution from rf^. 



is concluded to be zero. Thus the EAR spectrum in the 
FM phase shows no signal in the entire w regime. 

LP phase. — The MF spinor ^lp perpendicular to the 
magnetic field leads to Dk = aok, S% — (aifc + a_ifc)/-\/2, 

^k, ~ i{o-ik — a-ifc)/>/2, and 5^ = 0. The diagonalized 
Bogoliubov Hamiltonian is given [l^ as 



H. 



LP 
cff 



E 



Ekd\k)d{k) 



Eifl{k)Uk) 



(12) 



where Ei- 



\/£k{ek + 2nco) and 



\J i^k + '?)(efc + 9 + 2nci) are, respectively, a gapless- 
phonon and doubly-degenerate gapful-spin mode. The 
Bogoliubov transform gives the following form of the QZ 
perturbation as 



k^O i^=x,y 



/t(fc)/t(-fc) + H.c. 



(13) 



where -^q contains the terms commutable 
with the Hamiltonian (HI]), and = 

\J{ek + q + nci)2 - [eIY/2EI. The QZ perturba- 
tion accesses the two spin modes. From Eq. (IT3t . 
the retarded correlation function is straightfor- 
wardly calculated, and consequently the EAR spec- 
trum per unit volume is analytically obtained as 
i?LP(a;) = {6qfn^eYi{\Co\ - 2A5~)rLP(cj) with 



rLP(w) = 2\ 



2g- sgn(ci) 



(14) 



where fe(x) is a Heaviside's step function, and TZc — 
ri(2M|ci|)'^/^/647r?i^ with the system volume is a con- 
stant. Here, the following normalizations have been used: 
Cb — huj/2n\ci\, q = g/2n|ci|, and A, = ^q{q + sgn(ci)) 
which means a spin gap normalized by 2nci. Equa- 
tion p4p for various values of q is plotted in Fig.[21[a). We 



note that the gap of the EAR spectrum, 2Aq, closes on 
the phase boundaries, q = Q with ci > and q = 2n|ci| 
with ci < 0. 

TP phase. — The MF spinor ^tp provide the fluctu- 
ations Dk = [aik + a_ifc)/V2, = aok, = 0, and 
5"^ — (aifc — a_ifc)/\/2, and the Bogoliubov Hamiltonian 
is diagonalized 16] as 



ttTP 
^eff 



J2 [EkdHk)d{k) + EUl{k)h{k) + Elf\k)f{k) 



k^O 



(15) 



where Ek and E^ = ek{(-k + 2ci) are gapless modes of 



phonon and spin, respectively, and E^, is another spin 
mode with the gap A^. The QZ operator is written as 



X-TP 



_ -FTP 



k^O 



fik)f\-k) + R.c. 



(16) 



where -^q^ 
with the 



commutable 

Ar = 



denotes the terms 
Hamiltonian (jisp , and 

{ek + \q\ + ci)^ - {El)y2El. From Eq. 1^, the 
QZ modulation is found to access only the gapful-spin 
mode. Straightforwardly the retarded correlation func- 
tion of Eq. ([T5| is calculated and the EAR is analytically 
obtained as i?'^P(w) = ((5g)27^c6lH(|w 



-2Ag)rTP(w) with 



„TP 



(^) 



-m 



1 



(17) 



which is illustrated in Fig. [51 

It is remarkable that in the positive ci regime we have 
i?^P(w) = 2i?'rP(w) for the fixed \q\, and the factor 2 
comes from the accessible number of the spin modes by 
the QZ perturbation; namely the perturbation ([T5)) for 
the TP phase involves the two gapful-spin modes, while 
only one gapful-spin mode is accessed for the LP phase. 
Thus, the qualitative difference of the spectrum intensity 
allows us to distinguish the two polar phase. In addition, 
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for this different spectrum intensity we find the following 
result: if we continuously change q across q — 0, the dis- 
continuous spectrum change is observed. This would be 
interpreted to be a signal of the first order phase transi- 
tion between the polar phases. 

In summary, the EAR in the polar phases has two im- 
portant features: The first is that we can measure the 
spin gap Aq in the polar phases. In particular this spin 
gap is indeed a dominant feature of the low-energy spin 
excitation. The second is that by the discontinuous dif- 
ference of the spectrum intensity, we can observe the first 
order phase transition from the dynamical viewpoint. As 
we will discuss later, these conclusions do not change even 
if we have trap potentials. 



BA phase. — The MF spinor ^ba provides Dk 
sm9{aik + a-ik)/V2 + cosOaok, ~ cos9{aik 

a_ifc)/V2 -I- sin6'aofc, 5*1 = icos6{aik ■ 



^ik)/V2, and 



= sin0(aife — a_ifc), and under the condition Cp ^ |ci | 
the Bogoliubov Hamiltonian is diagonalized [iT 
as 



ttBA 
^cff 



E 



E^^ {k)dHk)d{k)+Ef^{k)fl{k)Mk) 



E%^{k)fly{k)Uy{k) 



(18) 



where Ef'^ik) y /efcjefc + 2n(co - Q|ci|)], 

E f^jk) ^ Vek{ek + q) and E^^{k) 

(efc -I- 2n|ci|)[efc -I- 2(1 — g2)n|ci|] are interpreted 
to be the density, the gapless- and gapful-spin mode, 
respectively. In addition, the QZ perturbation is 
represented as 



+ cT!LmL{-^) + vl^]Uk)d\-k) + H.c. 



where T, 



BA 



(19) 



commutes with the Hamilto- 



nian (TTSt . and the factors are 
x/(efc + ~qn\c^\Y - (£;?A(fc))2/2£;BA(fc)^ 



BA 



sin 26 



4 

sin 29 



•)BA 



{s/okPk 



] with ak = E^Mk)l{ek + 2n|ci|) 



4 LV"«f^K ^^^^i 

and Pk = E^^'^'-'^'^/ek- From Eq. (HH), the EAR spec- 
trum is given as R^^{u}) = {6q)^nc[rf^iu}) + 6'H(|tD| - 



2ABA)rBA(^) + 0^(\C,\ - AfA)^BA(^)J ^ith 



„BA 



„BA 



n2 



Cb'^ ^ q 



q^ +0)2 



rr(-) = (2|£i|)^/^ (Af)'l-I 



X (PI - Af ) 



where 7(0;) 



:|i;|-A^'A)[|5,|2(|^|_ABA)2 + |g^|] 



(20) 



(21) 



, (22) 



A?A ^ 



-\/l — g2 denotes the normalized gap of the spin mode 
E^^, and the normalized coupling ci = ci/cq has been 
introduced. 

Equation (|22l) for the various g's is shown in Fig. [2Jc) . 
The gapless spectral weight rf^{u}) describes a two- 
particle excitation of the gapless-spin mode Ef-^{k), and 
at the limit g 1 it is identical to R'^^{uj) at g — >■ 1. The 
weight r^y{i^), which is the two-particle excitation of the 
spin mode Ef^{k), gives the gapful EAR spectrum with 
the gap, 2A?A^ and quadratically vanishes at the limit 
g — > 0. The other spectrum weight r^A(^) from the 
pair excitation of the quasiparticles of the gapless-phonon 
Ef-^{k) and of the gapful-spin mode E^^{k) provides a 
gapful spectrum with the gap A?A, xhe two-particle 
excitation of the spin and density mode is peculiar to 
the BA phase as seen in the form of the QZ perturba- 
tion (IT9I) . In addition, the spectrum weight r^-^{ijj) is 
of order a/cq/ci while the others are independent of cq. 
Thus the EAR for Cq ^ |ci| is dominated by r^^{uj) in 
the frequency regime where r^^{ijj) is finite. 

Trap effect. — Based on the above results in the homo- 
geneous case, we discuss the EAR spectrum for trapped 
systems by local density approximation. Then, the spec- 
trum is calculated by taking the average of the local spec- 
trum with the weight of the local density n{r). The local 
EAR is obtained by replacing the mean density n by the 
local one n(r). Since the density n accompanies the in- 
teraction constant Cq and Ci in all the results, it turns out 
that the inhomogeneity modifies the strength of the inter- 
actions: The effective interaction around the trap center 
behaves stronger, and weaker as going away from the cen- 
ter. Thus, if the trap center is in the polar phase and in 
the FM, the whole system exhibits polar and FM state, 
respectively, as expected from Fig. [TJ In addition, since 
the gap of the EAR spectrum is independent of the den- 
sity and the interaction, the gapful feature of the EAR 
in the polar phases should remain even in the trapped 
case. On the other hand, when the trap center is in the 
BA phase, the outward regime would be polar state from 
the phase diagram [TJ However, since the EAR in the 
BA phase is characterized by the gapless spectrum, the 
qualitative feature is expected to be protected. From the 
above discussion, the EAR spectrum can be concluded 
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to characterize the phases of spin-1 BECs regardless of 
the presence of trap potentials. 

Conclusion. — We have discussed the response of the 
spinor BECs to modulation of the magnetic field by using 
linear response theory, and found the EAR spectrum to 
access the correlation of the QZ term. By calculating 
the EAR of the spin-1 BECs, we have demonstrated the 
spectrum as a function of the modulation frequency to 
show the explicit difference in each phase. In addition, we 
have pointed out that our results are robust against the 



trap effect. We note that while we have considered only 
spin-1 BECs in this paper, application to spin-2 BECs is 
straightforward. 
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